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Abstract 

Transport aircraft performance is strongly influ- 
enced by the effectiveness of high-lift systems. De- 
veloping wakes generated by the airfoil elements 
are subjected to strong pressure gradients and can 
thicken very rapidly, limiting maximum lift. This 
paper focuses on the effects of various pressure gra- 
dients on developing symmetric wakes and on the 
ability of a linear eddy viscosity model and a non- 
linear explicit algebraic stress model to accurately 
predict their downstream evolution. In order to re- 
duce the uncertainties arising from numerical issues 
when assessing the performance of turbulence mod- 
els, three different numerical codes with the same 
turbulence models are used. Results are compared 
to available experimental data to assess the accuracy 
of the computational results. 

1 Introduction 

Within the airframe industry, the optimal design 
of efficient high-lift devices for take-off and land- 
ing conditions is an important issue and one of the 
most challenging problems. For the current gener- 
ation of commercial aircraft, such high-lift systems 
require a simplicity of design but a high level of effi- 
ciency. Achieving this goal demands improved Com- 
putational Fluid Dynamics (CFD) design tools that 
can realistically model the high-lift system flow field 
for full scale free flight Reynolds numbers. Diffi- 
culties are inherent to both the geometric complex- 
ity of multi-element airfoils as well as limitations 
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in flow physics modeling. The flow field surround- 
ing a multi-element airfoil is dominated by very 
complex viscous flow phenomena such as boundary 
layer transition, laminar separation bubbles, conflu- 
ent turbulent boundary layers, multiple viscous wake 
interactions, and possible shock/boundary layer in- 
teractions. For this reason, the demands are high on 
the numerical methods used in this type of calcula- 
tion, and several crucial aspects must be treated in 
an adequate way. One of these aspects is turbulence 
modeling. 

In recent studies, 1,2 analyses of high-lift multi- 
element airfoil configurations were performed to as- 
sess the predictive capability of different types of 
turbulence models. The studies revealed three ar- 
eas where turbulence model predictions were possi- 
bly deficient and which affected the overall predic- 
tion of the flow field. These areas were the inability 
to predict transition location, the inability to pre- 
dict the slat wake, and the inability to properly ac- 
count for streamwise curvature effects in the main 
element /flap region. Using the multi-element airfoil 
configuration, it is difficult to separate these effects. 
In this study, the dynamics of a developing symmet- 
ric wake generated by a splitter plate is studied by 
focusing on the ability of two turbulence models to 
predict its downstream evolution in different types 
of pressure fields. While the majority of previous 
studies have focused on fully developed wakes far 
downstream of their origin, 3 the emphasis here is 
the “near” wake region downstream of the slat. 

Experimental studies which have isolated pres- 
sure gradient effects on a developing wake are rather 
scarce. 4 In this paper, numerical results are com- 
pared to the experimental results obtained at the 
University of Notre Dame (USA). 5 This experiment, 
which isolated the effect of pressure gradient on the 
development of a 2-D wake, has been specially de- 
signed to provide high quality data for use by the 
CFD community to validate turbulence models for 
wake flows. 

In addition to having a well-documented experi- 
mental study, it is of equal importance to be able to 
evaluate the turbulence models within a numerical 
code whose numerical characteristics are well known. 
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Even then it is optimal, although difficult, to per- 
form the validation studies using different numeri- 
cal codes in order to quantify any numerical biases. 
The challenge in this regard is then to be sure that 
each turbulence model is similarly implemented into 
the different codes. This level of effort is time con- 
suming, but becoming more essential in the current 
technical environment where an increasing number 
of ill-posed validation studies are being undertaken. 

The purpose of this paper is to describe a valida- 
tion study in which three different numerical codes 
CFL3D, 6 NSMB, 7 and PAB3D 8,9 are used to vali- 
date two different types of turbulence models on a 
near wake flow in zero, adverse, and favorable pres- 
sure gradients. The two different types of turbu- 
lence models include a linear isotropic eddy viscosity 
model 10 and an explicit algebraic stress model. 11 

The next section describes briefly the experimen- 
tal study performed at the University of Notre 
Dame. Then, the turbulence models are presented 
and the main differences between the three codes 
are highlighted. For each code, the linear isotropic 
eddy viscosity model and the explicit algebraic stress 
model are first validated against experimental data 
obtained in a simple channel. The last part of the 
paper presents the results obtained by the three 
codes using these two types of turbulence models 
for wake flows subjected to different types of pres- 
sure gradients. 

2 Experimental Data 

The data used in this study is the result of an 
experimental research program on wake develop- 
ment in different pressure gradients. This work is 
a collaborative effort between the Hessert Center for 
Aerospace Research at the University of Notre Dame 
(USA) and the NASA Langley Research Center 
(USA). With this carefully controlled experiment, 
it has been possible to isolate the effect of pressure 
gradient on the development of a 2-D wake in condi- 
tions similar to that experienced by the slat wake of 
a multi-element airfoil. One of the unique features of 
this experiment was the elimination of unnecessary 
complicating factors such as wake/boundary layer 
interaction and cove separation, which are present 
on real airfoil and wing configurations and can make 
meaningful CFD assessments an extremely difficult 
task. 

The two-dimensional wake experiments were per- 
formed in a low-speed wind tunnel facility. The 
Reynolds number was approximately 2.4 x 10 6 based 
on the length of the wake-generating body. For com- 
parison, it may be noted that a Boeing 737-100 op- 
erating at a Reynolds number (based on mean av- 


erage chord) of 15.7 x 10 6 during landing approach 
has a Reynolds number of about 1.8 x 10 6 based 
on slat chord. 5 The wake-generating body is a two- 
dimensional splitter plate with round nose and ta- 
pered trailing edge. The boundary layer on the flat 
plate is tripped by distributed roughness over the 
nose. The splitter plate (1.2192 m long) is used to 
generate the wake. The wind tunnel section with 
adjustable tunnel walls located immediately down- 
stream of the wake-generating body is used to pro- 
duce the desired pressure gradient environment for 
wake development. By using the adjustable walls, 
the initial degree of asymmetry can also be con- 
trolled. To date, the experimental study has been 
limited to the study of the structure and develop- 
ment of symmetric wakes. 

For different combinations of imposed pressure 
gradients, experimental mean velocity and turbu- 
lence data were compiled using hot wire and two- 
component Laser Doppler Velocimetry measure- 
ments. Experimental results have verified the two- 
dimensionality of the flow Held. They have also 
shown the effects of favorable and adverse pressure 
gradients on wake thickening and velocity defect de- 
cay. Note that a Reynolds stress bias error discov- 
ered in the u'v' measurements in Liu et al 5 has been 
corrected. See Liu 12 for more details. 

3 Numerical Approach 

3.1 Turbulence models 

Two-equation turbulence models represent the 
simplest level of closure that contains transport 
equations for both the characteristic velocity and 
length scale. However, it is well known that mod- 
els based on the Boussinesq eddy viscosity concept 
lack the universality required as they provide a some- 
what restrictive representation of the physical pro- 
cesses governing the Reynolds stress tensor. Typi- 
cally, they do not adequately account for the stress 
anisotropies induced by the flow. The Differen- 
tial Stress Model (DSM) devised by reference to 
the exact but unclosed equation for the Reynolds 
stress tensor can properly account for the advective 
and diffusive transport of the Reynolds stress ten- 
sor but is much more complicated and too expen- 
sive for many practical engineering applications. As 
an alternative, explicit algebraic stress models have 
been developed. They provide a close coupling with 
the DSM in that they are directly derivable from 
the stress transport equations subject to constraints 
on the advection and diffusion terms. Thus, these 
models retain the characteristics of the full DSM 
introduced through the pressure-strain correlation 
model but also bring some economy to the compu- 
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tation. In this study, two turbulence models have 
been selected: a linear, isotropic eddy-viscosity K- 
s model (EVM), 10 and an explicit algebraic stress 
model (EASM). 11 

3.2 Linear Eddy Viscosity K-e model 

The linear eddy viscosity model used in this study 
is based on the K — £ model of Speziale and Abid. 10 
The modeled transport equations are the following: 


where S (= Sij) and W (= Wij) are the strain rate 
and rotation rate tensors, respectively, 


, _ 1 f dui_ 

' ij ~ 2 \dxj 


du£\ = 1 fdui_ _ duj \ 

dxi ) ’ u 2 \ dxj dxi J 


and the oq are scalar coefficient functions of the 
invariants rf (= {S 2 }) and 1Z 2 (= — {W 2 }/{S 2 }). 
The coefficients oq, a 2 , and 03 are given by 11 
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where D / Dt = d/dt + Uid/dxi, V is the produc- 
tion of turbulent kinetic energy (= —Tijdui/dxj), 
Tij is the Reynolds stress tensor (= —v t {dui/dxj + 
duj/dxi), K(= 7jT nn ) is the turbulent kinetic en- 
ergy, e is the turbulent dissipation rate, and v t is 
defined as v t = C^K 2 /e, where = 0.09. 

The role of the function f E is to remove the singu- 
larity in the destruction of the dissipation term; its 
expression is given by: 
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where ao = (l — | a\a\ri 2 + 2a2a|?y 2 7?. 2 ) and r = 
K/e. When R = cqS, then {RWS} = {RS 2 } = 0 
and {RS} = airf 2 . The coefficients a; are directly 
related to the pressure-strain correlation model used 
in closing the stress transport equation. This study 
uses the SSG pressure-strain model, 16 which yields 




02 = -(2-C-4) . 


fe = l-exp(-Rex/12), (3) 

with Rex = s/Kd/v the Reynolds number based on 
the distance to the nearest wall. 


3.3 Explicit Algebraic Stress Model 

The development of the explicit algebraic stress 
model (EASM) used in this study is outlined in de- 
tail in Rumsey et al. 11 It is based on the formu- 
lation developed by Jongen and Gatski 13 from the 
earlier work of Pope 14 and Gatski and Speziale. 15 
The reader is referred to these earlier studies for ad- 
ditional background. The resulting EASM for the 
anisotropy tensor 
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in two-dimensional flows is given by 
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Cl = 3.4( 7l = 0.7), Cl = 1.8(70 = 1.9), 

C' 2 = 0.36, C 3 = 1.25, C 4 = 0.4. 

The explicit tensor representation given in Eq. (5) 
is coupled with the K-e two-equation model de- 
scribed in Eqs. (l)-(3). Only the coefficients C M and 
f E require slightly different values and are now given 

by 


b = aiS+a 2 (SW-WS)+a 3 
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^{S 2 }l) , (5) 


CV = 0.096, 


f E = 1 - exp 
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In the development of the algebraic stress model, 
two assumptions are made in the formulation of the 
implicit algebraic equation for the stress anisotropy. 
The first is an equilibrium assumption on the Rey- 
nolds stress anisotropy 

Dbjj _ Drjj _ Tjj DK 
Dt Dt K Dt 

and the second is an assumption on the anisotropic 
behavior of the turbulent transport and viscous dif- 
fusion term T>ij , that is 

Vij - = 0. (14) 

Equation (14) assumes that the anisotropy in T>ij 
is directly related to the anisotropy in the Reynolds 
stresses themselves. However, Gatski and Rumsey 17 
has shown that near the center region of wake flows 
Eq. (14) leads to unphysically large levels of turbu- 
lent eddy viscosity. A different constraint on T>ij 
can be found by rewriting Eq. (14) in terms of the 
anisotropy tensor bij. This leads to the modified 
constraint 

T>ij - \v nn 8ij = 0, (15) 

which simply states that the deviatoric part of the 
tensor T>ij is zero. The contribution from bijV nn 
now appears as a modification to the a 4 coefficient 
in the form 
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As is well known, the algebraic stress formula- 
tion should yield the same results as a full second- 
moment closure in the case of homogeneous shear. 
Thus a modification to the algebraic stress model 
as described here needs to maintain this relation- 
ship. This requirement can be fulfilled by invok- 
ing an equilibrium assumption on the turbulent time 
scale, that is 


Dt 1 DK K De 
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For the case of homogeneous shear, this yields the 
relationship 
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which can be substituted into Eq. (17) along with 
the equilibrium condition on V /£, 
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The values for C s \ and C s 2 are the same as those 
used in the linear EVM, that is, 1.44 and 1.83, re- 
spectively. This modification meets the requirement 
that the homogeneous shear results are unaltered, 
and only slightly affects the results for the log-layer 
where the value of C M is changed from 0.096 (see 
Eq. (12)) to 0.0885. The only other alterations to 
the algebraic stress formulation are that now a\, Jq, 
and 7* are used instead of the coefficients 04, 70, 
and 71 . 


3.4 Numerical Algorithms 

In order to provide a controlled assessment of 
the performance of the turbulence models, they 
have been implemented into three different numer- 
ical codes CFL3D, 6 NSMB, 7 and PAB3D. 8 - 9 With 
this multiple-code implementation, the effect of the 
purely numerical influences on the solution results 
can be quantified. Many organized efforts at val- 
idating the performance of turbulence models suf- 
fer from the lack of precise coordination of detailed 
turbulence model implementation. In the valida- 
tion study presented here, this deficiency is recti- 
fied and a more accurate assessment of the turbu- 
lence models can be made. The three codes solve 
the three-dimensional Reynolds-averaged compress- 
ible Navier-Stokes equations. The distinctive fea- 
tures of each code are briefly described in this sec- 
tion. 


CFL3D 

CFL3D is a computational method for structured 
grids developed at NASA Langley Research Cen- 
ter. The code uses an upwind Unite volume for- 
mulation. It neglects viscous cross-derivative terms, 
which results in the thin-layer Navier-Stokes equa- 
tions in specified coordinate directions. Third-order 
upwind-biased spatial differencing on the convective 
and pressure terms, and second-order differencing 
on the viscous terms are used; it is globally second- 
order spatially accurate. The CFL3D code can solve 
flow over multiple-zone grids that are connected in 
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a one-to-one, patched, or overset manner, and can 
employ grid sequencing, multi-grid, and local time 
stepping when accelerating convergence to steady 
state. Upwind-biased spatial differencing is used for 
the inviscid terms, and flux limiting is used to ob- 
tain smooth solutions in the vicinity of shock waves, 
when present. No limiter was employed for this 
study. Viscous terms are centrally differenced. The 
flux difference-splitting (FDS) method of Roe is em- 
ployed to obtain fluxes at the cell faces. The CFL3D 
code is advanced in time with an implicit three- 
factor approximate factorization method. The im- 
plicit derivatives are written as spatially first-order 
accurate, which results in block-tridiagonal inver- 
sions for each sweep. Flowever, for solutions that 
utilize FDS the block-tridiagonal inversions are fur- 
ther simplified with a diagonal algorithm. Turbu- 
lence equations are solved uncoupled from the mean 
equations. 

NSMB 

The NSMB (Navier-Stokes Multi-Block) code is 
a finite volume code for structured grids devel- 
oped by two universities, EPFL (Switzerland) and 
KTF1 (Sweden), a research institute CERFACS 
(France), and two companies, Aerospatiale (France) 
and SAAB (Sweden). A cell-centered finite vol- 
ume method is used, and the viscous cross-derivative 
terms are included in the formulation. The equa- 
tions are discretized in space using either a central 
difference scheme yielding second order accuracy in 
space for a smooth grid with artificial dissipation 
(standard Jameson dissipation, Martinelli dissipa- 
tion or matrix dissipation), or an upwind scheme 
(Roe and AUSM using different limiters). In this 
study, the central scheme with the standard Jameson 
dissipation has been used. The time integration can 
be achieved through the use of an explicit Runge- 
Kutta time stepping. For steady-state problems, an 
implicit time integration scheme is available in the 
form of the LU-SGS scheme of Yoon and Jameson 18 
which is based on a Lower-Upper factorization and a 
symmetric Gauss Seidel relaxation. The turbulence 
equations are solved coupled to the mean equations. 

PAB3D 

PAB3D 8,9 is a code developed by Analytical Ser- 
vices and Materials, Inc. for structured grids us- 
ing an upwind finite volume formulation. It can ei- 
ther solve thin-layer Navier-Stokes or any two grid- 
index directions with fully coupled viscous terms. In 
this study, the two directions were uncoupled as in 
CFL3D. The Roe FDS scheme is used to evaluate the 
explicit part of the governing equations and the van 


Leer scheme is used to construct the implicit opera- 
tor. The FDS solutions with this method were not 
further simplified with a diagonal algorithm which 
should only impact solution time and not solution 
accuracy. The diffusion terms are centrally differ- 
enced, and the inviscid flux terms are upwind dif- 
ferenced. Similar to CFL3D, several different flux 
limiters can be applied. Current solutions used the 
min-mod flux limiter. Both the mean and turbulent 
flow equations used a third-order-accurate scheme 
for the convective terms and second-order accurate 
for the viscous diffusion terms. The turbulent flow 
equations are solved uncoupled from the mean flow 
equations. 

4 Results and Discussion 

4.1 Turbulence model validation 

As a first step in the process of turbulence model 
calibration and validation, all three codes CFL3D, 
PAB3D, and NSMB used the two-equation turbu- 
lence model (EVM) and explicit algebraic stress 
model (EASM) described in the previous section to 
compute the channel flow of Laufer. 19 The same grid 
was used in all three codes. The grid models half the 
channel and consists of 89 x 73 gridpoints. The cells 
are clustered close to the wall in the y direction and 
equally spaced in the x direction. The mean value 
of y + at the first gridpoint off the wall is below 1. 

The free-stream Mach number is set to 0.2, and 
the Reynolds number based on the freestream veloc- 
ity and half of the channel width is 61600. The total 
pressure and total temperature are prescribed at the 
inlet. A back static pressure is imposed at the outlet. 
As shown in Figs. 1 and 2, both mean flow and tur- 
bulence quantities are in a very close agreement for 
the three CFD codes when either the EVM or EASM 
are employed. In general, both turbulence models 
yield results in good agreement with experiment. As 
expected, the eddy viscosity model cannot predict a 
difference between the normal Reynolds stress com- 
ponents (Fig. Id), whereas the normal stress com- 
ponents predicted by the EASM (Fig. 2d) show the 
correct trends. 

4.2 Wake Flow 

4.2.1 Numerical implementation 

The computational domain for the wake flow test 
case (zero pressure gradient) is depicted in Fig. 3. 
The domains for the adverse and favorable condi- 
tions have different tunnel wall shapes downstream 
of the splitter plate. The overall geometry is iden- 
tical to the experimental one, except that an ellip- 
tic nose is used for the splitter plate instead of a 
circular nose. Initial computations with a circular 
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nose showed a large laminar separation bubble that 
was quite difficult to compute and some unsteadi- 
ness was noticed. Because the experimental tests 
included roughness over the nose to trip the bound- 
ary layer on the plate, the laminar bubble probably 
did not exist in the experiments, although the exper- 
imental data were insufficient in this region to verify 
this. The use of an elliptic nose in the computational 
study eliminates the laminar bubble. 

Consistent with experimental conditions, no-slip 
and adiabatic wall conditions are assumed at solid 
walls. At the inflow boundary, the total pressure 
and total temperature are specified, and at the out- 
flow boundary, the back static pressure is specified. 
As in the experimental study, zero, adverse and fa- 
vorable pressure gradients are attained in the nu- 
merical study by using a different geometry for the 
tunnel wall shape. It was necessary to test different 
back pressures for each of the test cases to recover 
the experimental velocity at the end of the splitter 
plate; these tests have been carried out essentially 
with CFL3D. Then, for each test case, the ratio of 
the static pressure at the outflow P ou tfiow over the 
total pressure at the inflow Pt, nftow has been fixed at 
the same value in all three codes. Because PAB3D is 
not a non-dimensional code and because computa- 
tions for these test cases with NSMB used physical 
values as well, the pressure and temperature ratios 
used in CFL3D were converted to equivalent met- 
ric quantities assuming a freestream total pressure 
of 101.325 KPa and total temperature of 300 K. 

Because the wakes are assumed symmetric, 5 com- 
putations are only performed on the upper half of 
the computational domain and a symmetric bound- 
ary condition is imposed in the center of the tunnel. 
This minimizes grid points and allows for a greater 
concentration of the points in regions of interest. An 
early study with CFL3D confirmed that identical re- 
sults were obtained using a half domain or full one. 

Computational grids, comprised of 5 blocks, used 
a high concentration of points near the walls, plate 
surface and in the wake. All grids were generated 
with an average distance from the wall to the first 
gridpoint of approximately y + = 1.1. The total 
number of points in the finest mesh is 63549. 

Grid dependence tests have been performed using 
CFL3D on three different grids (each successively 
coarser grid made up of every other point from the 
finer levels). Of the cases tested, the favorable pres- 
sure gradient test case showed the greatest sensitiv- 
ity to grid size. The velocity and turbulent shear 
stress sensitivities for this case using the EVM can 
be seen in Figs. 4 and 5. The trends suggest that 
further mesh refinement may be useful, but should 


not yield significantly different results from the fine 
grid. Although not shown, results for the zero pres- 
sure gradient test case showed almost no sensitivity 
in the velocity profiles. Sensitivities using EASM 
were similar to EVM. 

Because the turbulence models are not able to 
accurately predict transition location, the point of 
transition was different for each model as well as 
for each code. In all the tests, the codes were run 
“fully turbulent” : transition location was not fixed 
but the respective models and codes were allowed to 
transition freely. Resulting transition locations for 
both EVM and EASM are summarized in Tables 1 
and 2, respectively, for the adverse pressure gradient 
(APG), zero pressure gradient (ZPG), and favorable 
pressure gradient (FPG) cases. In general, PAB3D 
tends to transition somewhat further downstream 
than CFL3D; however, as will be shown later, these 
differences do not lead to significant differences in 
the computed solutions at the end of the plate. With 
the exception of FPG using EVM, NSMB’s transi- 
tion locations are consistently further downstream 
than the other codes. 



CFL3D 

NSMB 

PAB3D 

APG 

0.053 

0.370 

0.368 

ZPG 

0.064 

0.355 

0.206 

FPG 

0.150 

0.057 

0.226 


Table 1. Transition ^-location in meters along the 
splitter plate using EVM. 



CFL3D 

NSMB 

PAB3D 

APG 

0.056 

0.370 

0.333 

ZPG 

0.099 

0.345 

0.213 

FPG 

0.191 

0.310 

0.180 


Table 2. Transition ^-location in meters along the 
splitter plate using EASM. 


4.2.2 Zero Pressure Gradient 

The zero pressure gradient wake test case is com- 
puted to serve both as a baseline for comparison with 
variable pressure gradient wake development and as 
a comparative check on the predictions of the three 
codes using the same turbulence models. 

The symmetric turbulent wake generated by a 
splitter plate is typically divided into three regions. 
The first region begins at the end of the plate and 
may extend to approximately 25 momentum thick- 
nesses, which corresponds in this study to the lo- 
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cation x = 1.3 m. The second, or intermediate, 
region modifies the logarithmic region of the plate 
boundary layer. The third, or far wake, region be- 
gins about 350 momentum thicknesses downstream 
of the trailing edge where the mixing is complete. 20 
This last region is located at x = 2.3 m in this study. 

As an initial comparison of code and model be- 
havior, predictions from the three codes were com- 
pared to experimental mean velocity and shear stress 
data at x = 1.2167 m, near the end of the splitter 
plate. The velocity profiles (see Fig. 6) predicted 
by CFL3D and PAB3D show good agreement with 
experimental data for both the EVM and EASM, 
whereas NSMB leads to a fuller profile and a smaller 
boundary layer thickness in both cases. The shear 
stress profiles at the same location (see Fig. 7) show 
greater variation in the case of EASM, and the agree- 
ment with experiment is not as good as the mean 
velocities for either turbulence model. 

Because great pains were taken to ensure the same 
implementation of the turbulence models in all three 
codes (recall the channel flow validation study), the 
variations in the results for this more complex case 
are somewhat disappointing. One possible cause 
for the variations may be the different numerics in 
the three codes: NSMB uses a central scheme with 
the standard Jameson dissipation and CFL3D and 
PAB3D use an upwind scheme. However, while 
a grid study with NSMB (not shown) indicates a 
somewhat greater sensitivity to grid than CFL3D, 
its trend shows decreasing boundary layer thickness 
with increasing grid density. This trend does not 
account for the differences seen in Fig. 6. 

Another possible cause for the different results 
is the fact that the three codes transition in dif- 
ferent locations. In Tables 1 and 2, the transition 
locations for ZPG are successively further down- 
stream (by between 0.1 and 0.15 m) for the three 
codes. This difference does not seem to affect the 
results between CFL3D and PAB3D, but it is pos- 
sible that NSMB’s transition location is far enough 
downstream to cause a noticeable difference in re- 
sults near the trailing edge. 

The computed mean velocity profiles in the wake 
are compared with the experimental results in 
Figs. 8 and 9. In the near wake (Figs. 8a and 9a), the 
agreement between numerical and experimental ve- 
locity profiles is still very good overall; whereas, fur- 
ther downstream (Figs. 8b, c and 9b, c) the velocity in 
the wake is slightly overpredicted by all codes. The 
initial fuller profile displayed by the NSMB code per- 
sists, but the qualitative trend is the same as both 
the CFL3D and PAB3D codes. 


Figures 10 and 11 show the comparison between 
the computed shear stress and experimental data. In 
the near wake region (Figs. 10a and 11a), the com- 
puted results are in good overall agreement with the 
experimental data. Further downstream in the wake, 
the computed shear stresses (Figs. 10b, c and llb,c) 
decrease too fast compared to the experimental data. 
It is interesting to note that the qualitative trends 
are not dependent on the choice of turbulence model. 
Thus, the added physics associated with the EASM, 
such as better accounting for stress anisotropy and 
closer sensitivity to gradients of the mean velocity, 
do not play a significant role in the prediction of the 
developing wake flows. 

It should be noted that while transition location 
and numerical differences (such as central differenc- 
ing vs. upwind differencing) are likely to be the cause 
of most of the differences observed between the three 
codes, various implementation details in the codes 
are also very different and may be factors. Some 
examples include the order of accuracy of the tur- 
bulent advective terms, type of flux limiter, method 
for including the nonlinear turbulent stress terms, 
and nondimensionalization relative to the boundary 
conditions. These types of differences were too nu- 
merous to explore for this study. 

4.2.3 Adverse Pressure Gradient 

The experimental results have shown that the in- 
fluence of the pressure gradient on wake development 
and structure is very significant. 5 When an adverse 
pressure gradient is imposed, the wake width is in- 
creased and the Reynolds shear stress is amplified 
(relative to ZPG). 

The nominal pressure gradient in the wake for this 
case is dC' p /dx = 0.3386/m. Numerical results for 
the mean velocity obtained with the three different 
codes are in good agreement with experimental re- 
sults for both the EVM (Fig. 12) and the EASM 
(Fig. 13) models. The computed turbulent shear 
stresses are shown in Figs. 14 and 15. The computed 
shear stresses over-predict the local maximum exper- 
imental data initially and then further downstream 
under-predict the experimental results. In compar- 
ison with the corresponding zero pressure gradient 
case (Figs. 8 - 11), the numerical predictions for the 
adverse pressure gradient case are better and the re- 
sults are more consistent between the different codes. 

4.2.4 Favorable Pressure Gradient 

When the wake develops in a favorable pres- 
sure gradient, the wake width is reduced and the 
Reynolds stress is decreased (relative to ZPG). Be- 
cause in all cases the numerical results for both the 
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EVM and the EASM are very similar, for brevity 
only the EASM results are shown here. 

The nominal pressure gradient in the wake for this 
case is dC p /dx = —0.5984/m. In Fig. 16 the devel- 
opment of the mean velocity profiles is shown. As 
in the adverse pressure gradient case, all three codes 
underpredict the wake thickness at x = 1.3462 m 
(Fig. 16a). The wake deficit is also slightly off. At 
the intermediate station (Fig. 16b) the computed re- 
sults continue to overpredict the mean velocity, while 
at the last station (Fig. 16c) the computed results 
are becoming more in line with the experiment. The 
computed turbulent shear stress profiles are shown 
in Fig. 17. The trends here are similar to the zero 
and adverse pressure gradient cases. 

5 Conclusions 

A purpose of this paper was to describe a study 
in which three different numerical codes, CFL3D, 
NSMB, and PAB3D, are used to validate two differ- 
ent types of turbulence models on a near wake flow in 
zero, adverse, and favorable pressure gradients. The 
two different types of turbulence models included a 
linear isotropic eddy viscosity model (EVM) and an 
explicit algebraic stress model (EASM). 

In order to perform as credible a validation study 
as possible, the turbulence models were implemented 
as identically as possible in each code. Then, for 
validation of each code, the EVM and EASM results 
were compared to experimental data obtained in a 
simple channel. Results showed excellent agreement 
among the three codes. 

The three codes were then applied to the study 
of symmetric turbulent wakes. Comparisons be- 
tween the numerical results and the experimental 
data yielded the following conclusions: 

• Mean velocity profiles for the three different 
pressure gradient conditions are, in general, in 
good agreement with the experimental profiles 

• More significant discrepancies are observed for 
the turbulent quantities, although both the 
EVM and EASM are able to qualitatively re- 
produce the experimental results 

• Although small differences exist, the EVM and 
EASM lead to very similar results 

The differences between computed results from 
the three codes for the wake cases cannot be at- 
tributed to the turbulence models, which are identi- 
cal in each code, but rather to the differences in tran- 
sition location and numerical differences between the 
codes. Further study into the influence of the numer- 
ical scheme as well as the influence of the artificial 


dissipation, especially on the turbulent stress results, 
is needed. A method for guaranteeing transition in 
a consistent fashion between different codes (when 
running “fully turbulent”) should be explored. 

As discussed in the Introduction, the motivation 
for this study was to try to isolate the predictive 
capability of turbulence models on developing wake 
flows which have relevance to the wakes shed from 
airfoil elements, such as slats, in high-lift systems. 
Such wakes develop in strong pressure gradients, and 
are an integral part of the aerodynamics of the high- 
lift system. The present study has shown that the in- 
clusion of the effects of turbulence anisotropy and/or 
improved sensitivity to the mean velocity gradient 
Held are not central to the accurate description of 
near wake development. Even a simple two-equation 
eddy viscosity turbulence model can reproduce the 
near-wake development fairly well in the three pres- 
sure fields studied. 

However, it should be noted that the pressure gra- 
dients experienced by a slat wake are typically vari- 
able (a short, very strong favorable gradient followed 
by an adverse gradient), and these gradients can 
be significantly stronger than those imposed in this 
wake study. A potentially useful future experimen- 
tal study would attempt to duplicate the variable 
pressure field levels seen by a typical slat wake. 

A final important conclusion from this study is 
that turbulence models can easily be (and often 
are) coded differently in different CFD codes. Small 
variations in the models (constants, damping func- 
tions, etc.) can completely alter computed results 
and muddy the conclusions. Ensuring identical tur- 
bulence model implementation is tedious; however, 
it is well worth the effort. Using multiple CFD 
codes for a turbulence model validation study ac- 
complishes the following. First, it gives confidence 
that each turbulence model has been implemented 
correctly, and, as a result, it is easier to draw conclu- 
sions regarding the turbulence model itself. Second, 
when transition location is the same, differences be- 
tween the codes’ results may be solely attributed to 
differences in the codes’ numerics; this is helpful for 
establishing limits on the numerical uncertainty for 
a given case. 
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Figure 1. Fully developed turbulent channel flow 19 
predictions with EVM: (a) Mean velocity, (b) 
Reynolds shear stress, (c) turbulent kinetic energy, 
(d) Reynolds normal stresses 


Figure 2. Fully developed turbulent channel flow 19 
predictions with EASM: (a) Mean velocity, (b) 
Reynolds shear stress, (c) turbulent kinetic energy, 
(d) Reynolds normal stresses 
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Figure 3. Computational domain for the zero pres- 
sure gradient wake test case 
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Figure 6. Wake flow with zero pressure gradient - 
Mean velocity on the splitter plate at x = 1.2167 m: 
(a)EVM, (b) EASM 


Figure 4. Effect of mesh refinement on the mean 
velocity for the favorable pressure gradient wake test 
case - EVM 
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Figure 5. Effect of mesh refinement on the turbulent 
shear stress for the favorable pressure gradient wake 
test case - EVM 


Figure 7. Wake flow with zero pressure gradient 
Turbulent shear stress on the splitter plate at x 
1.2167 m: (a) EVM, (b) EASM 
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Figure 8. Wake flow with zero pressure gradient - 
Mean velocity in the wake - EVM 



Figure 9. Wake flow with zero pressure gradient - 
Mean velocity in the wake - EASM 
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Figure 10. Wake flow with zero pressure gradient - Figure 11. Wake flow with zero pressure gradient - 
Turbulent shear stress in the wake - EVM Turbulent shear stress in the wake - EASM 
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Figure 12. Wake flow with adverse pressure gradient 
- Mean velocity in the wake - EVM 
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Figure 13. Wake flow with adverse pressure gradient 
- Mean velocity in the wake - EASM 
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Figure 14. Wake flow with adverse pressure gradient Figure 15. Wake flow with adverse pressure gradient 
- Turbulent shear stress in the wake - EVM - Turbulent shear stress in the wake - EASM 
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Figure 16. Wake flow with favorable pressure gradi- Figure 17. Wake flow with favorable pressure gradi- 
ent. - Mean velocity in the wake - EASM ent - Turbulent shear stress in the wake - EASM 
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